arXiv:1501.07856vl [astro-ph.HE] 26 Jan 2015 


Draft version February 2, 2015 

Preprint typeset using DTfTjX style emulateapj v. 2/16/10 


BUSTING UP BINARIES: ENCOUNTERS BETWEEN COMPACT BINARIES AND A SUPERMASSIVE 

BLACK HOLE 

Eric Addison 

Department of Physics, Utah State University, Logan, UT 84322, USA. 


Pablo Laguna 

Center for Relativistic Astrophysics and School of Physics, 

School of Physics, Georgia Institute of Technology, Atlanta, GA 30332, USA. 


Shane L. Larson 

Center for Interdisciplinary Exploration and Research in Astrophysics, 

Northwestern University, Evanston, IL 60208, USA 
and Department of Astronomy, Adler Planetarium, Chicago, IL 60605 USA. 

Draft version February 2, 2015 

ABSTRACT 

Given the stellar density near the galactic center, close encounters between compact object binaries 
and the supermassive black hole are a plausible occurrence. We present results from a numerical 
study of close to 13 million such encounters. Consistent with previous studies, we corroborate that, 
for binary systems tidally disrupted by the black hole, the component of the binary remaining bound 
to the hole has eccentricity ~ 0.97 and circularizes dramatically by the time it enters the classical LISA 
band. Our results also show that the population of surviving binaries merits attention. These binary 
systems experience perturbations to their internal orbital parameters with potentially interesting 
observational consequences. We investigated the regions of parameter space for survival and estimated 
the distribution of orbital parameters post-encounter. We found that surviving binaries harden and 
their eccentricity increases, thus accelerating their merger due gravitational radiation emission and 
increasing the predicted merger rates by up to 1%. 

1. INTRODUCTION 

Observations of gravitational waves (GWs) will allow 
us to probe dynamical astrophysical systems in regimes 
where strong gravity plays a key role. In galactic nuclei, 

GWs from stellar mass interactions with a supermassive 
black hole (SMBH) will be excellent probes of the prop¬ 
erties of the SMBH and the stellar population itself. It 
is thus important to investigate the types of encounters 
that should be expected as well as the strength and fre¬ 
quency of GWs emitted by these events. 

In the present work, we will focus on encounters 
between stellar mass binaries and a galactic SMBH. 

Previous studies have focused on the hyper-velocity 
star (HVS) produced by the d i sruption of the binary 
and stellar collisions dHills I 1 1988t i Cualandris et, al.ll2005l : 
lAntonini et all l2010f) . TMiller et al.1 (|2005li ~ investigated 
the eccentricity of the bound component created by the 
binary disruption, namely the extreme mass-ratio inspi¬ 
ral (EMRI) left behind. More recently, binary interac¬ 
tions with SMBH have been explored in the context of 
the GW emissions in the LIGO band from binaries driven 
to merger by Kozai resonance (jAntonini & Peretsll2012t 
lAntonini et al.|[2014f) . 

The strongest GW sources will involve compact, de¬ 
generate stellar remnants that can survive close encoun¬ 
ters with the central black hole, namely neutron stars 
and stellar mass black hole (BH)s. The fraction of field 
stars in binaries varies by stellar type, but for O- and 
B-type stars, which are the stellar types massive enough 
to form neutron stars or BHs, it is estimated that more 


than 75% of O-type an d 70% of B-type stars have some 
number of companions dRaghavan et al.ll2010D . Near the 
galactic center, the density of stars grows very large com¬ 
pared to field conditions, with densit y estimates up to 
10 8 Mq pc -3 within the inner O.lpc, (j Alexander!(2005 1. 
Given this information, it is reasonable to expect the ex¬ 
istence of compact object (CO) binaries in the galactic 
center, and in fact it has been observed by way of X- 
ray transients that CO binaries exist near the GC, and 
are m ore abundant within the inner 1 pc dMuno ct al.1 
12005 1. It has been estimated that as many as ~ 20,000 
stellar mass BH binaries have segregated withi n the inner 
» 1 p c of the Milky Way galactic center (lO’Learv ct al.l 
120091) . 

With the existence of CO binaries near the galac¬ 
tic center, it is reasonable to expect that some number 
of them may interact directly with the SMBH. Many 
known main sequence stars exi s t in bound o rbits around 
the SMBH (I Chez et al.l 120051 12008 1 120091) . suggesting 
that the same could be true for COs and CO binaries. 
Such interactions have implications for GW campaigns 
of observations, for example ground based observato¬ 
ries (e.g. LIGO) searching for compact binary coales¬ 
cences (CBCs). Proposed space based interferometers 
(e.g. LISA) will be sensitive to EMRIs which can be 
created by the tidal disruption of CO binaries near the 
SMBH. 

In the present study, we corroborate and extend the 
results of lMiller et al.l ([20051 1 for the formation of EMRIs 
by tidal disruption to arbitrary binary orientations. In 
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addition, we focus our attention on those binary systems 
that survive the encounter with the SMBH. In particular, 
we investigate whether those binaries are able to survive 
and face subsequent encounters before merging from GW 
emission or merge faster than those not experiencing en¬ 
counters with SMBH, and discuss the implication of the 
encounters on predicted CBC rates. We concentrate on 
initially circular, equal mass (mi = m 2 = m), CO bina¬ 
ries approaching the SMBH in a parabolic orbit. We vary 
the orientation of the orbital angular momentum of the 
CO binary relative to the angular momentum of the CO 
binary orbiting the SMBH. We also vary the pericenter 
distance r p between the CO binary and the SMBH. Our 
study explores a greater volume of the parameter space; 
lAntonini fe Per eta (|2012T ) focused o n binaries bound to 
the SMBH at a fixed orientation, and I Miller et all (|2005fl 
explored hyperbolic encounters with coplanar binaries. 
Considering a larger region of parameter space comes 
at the expense of integration sophistication, though the 
time and distance scales involved in our simulation sug¬ 
gest that relativistic effects will not play a significant role 
for the majority of the parameter range, and Newtonian 
gravitational forces should suffice. 

The paper is organized as follows: Section [5] outlines 
the outcomes of the three body encounters based on en¬ 
ergy arguments. In Section [3] we describe the setup of 
the encounters. Parameter probability distributions and 
an estimate of the tidal radius are given in section 0 
The analysis of the disrupted binaries is found in Section 
0 In Section 0 we discuss the results of those binaries 
that survived the encounter with the SMBH. Section [7] 
includes a discussion of the changes in the lifetime of 
the CO binary due to GW emission as a result of the 
encounter with the hole as well as consequences to GW 
detection rates for CO binaries. The paper ends with 
conclusions in Section 0 


proximate V 3 ft: 0. The total energy thus become 
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where V b is the velocity of the center-of-mass of the bi¬ 
nary relativity to the SMBH. Since for the incoming bi¬ 
nary R 3> r, with R denoting the distance of the center- 
of-mass of the binary to the SMBH, we can rewrite the 
total energy as 

I 9 GM h M. 

E ft -M b V b 2 -- + E b = E cm + E b . ( 6 ) 

z K 

As mentioned before, we inject the binaries in parabolic 
orbits; therefore, the initial center of mass energy of the 
binary with the black hole is E crn o = 0. Therefore, the 
total energy is given by the initial internal binding energy 
of the binary, E = E b o < 0. 

The effect of the encounter is to re-distribute the en¬ 
ergy available, i.e. E b fi, among the three bodies. If the 
binary survives, E b < 0. Given that 

E = E cm + E b = E b fi < 0 , (7) 

the CO binary after the encounter could be bound 
(E cm < 0) or unbound (E cm > 0) to the SMBH. 

On the other hand, if the CO binary does not survive 
( E b > 0 ), the separation r of its components will grow; 
thus, one can neglect in Eq. 0 
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2. ENERGETICS AT A GLANCE 

The total energy for a three-body system can be writ¬ 
ten as 
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In the barycentric frame, with M b = mi + m 2 the mass 
of the binary and M. = m 3 the mass of the SMBH, the 
total energy reads 
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where we have used again that r 3 = 0 and V 3 ft 0. There¬ 
fore, for CO binaries that are disrupted 

E = E\ + E 2 = E bt Q < 0 . (10) 

The possible outcomes in this situation are both COs 
bound to the SMBH (Ei,E 2 < 0), or one bound and 
the other unbound ( Ei < 0 < Ej). Clearly, the case in 
which bound components are unbound (0 < Ei,E 2 ) is 
not possible. 

We will refer to a CO binary in a bound orbit around a 
SMBH as a binary extreme mass-ratio inspiral (BEMRI), 
and will classify the aftermath of the binary from the 
encounter with the SMBH into one of the following four 
classes: 


Above, 
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is the internal energy of the binary with r = iq — r 2 , 
v = Vi — v 2 and /i = m\m 2 /M b . Furthermore, 


x 2 
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( 4 ) 


is the relative velocity between the binary and the 
SMBH. Since M. M bl we can set r 3 = 0 and ap¬ 


• DB: Disrupted binary. 

• Long BEMRI: Survived binary bound to SMBH 
with r gw > P,. 

• Short BEMRI: Survived binary bound to SMBH 
with r gw < P,. 

• SU: Survived binary unbound to the hole. 

In this classification, r gw refers to the binary merger life- 
time from G W emission, the so-called Peters lifetime 
(|Peterslll964f) . and P. is the period of the bound binary 
around the SMBH. One of the main motivations of our 
study is to investigate how the probability of a binary 
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falling into one of these categories could alter the pre¬ 
dicted CBC. 

The fate of a CO binary encountering a SMBH depends 
on its penetration factor (3. The penetration factor is 
defined as the ratio /? = r t /r p , where r t is the tidal radius 
and r p the distance of closest approach to the SMBH. In 
analogy with the common definition of tidal radius for 
the stellar disruption of stars by massive BHs, we define 
the tidal radius r* as 


n 





where ao is the initial semi-major axis of the CO bi¬ 
nary. The radius r t is an approximation to the distance 
where within the tidal forces by the SMBH exceed the 
self-binding energy of the binary. The exact distance will 
depend on the orbital parameter s and orientation of the 
binary (jHamilton fc Bumslll99lL 119921) . One of the ob¬ 
jectives of our study is to provide a statistical estimate 
of the tidal radius from our set of encounters. 

As mentioned before, the effect of the encounter is a 
redistribution of energy and angular momentum in the 
system. As a consequence, orbital parameters in the CO 
binary, such as eccentricity and semi-major axis are af - 
fected. Regarding eccentricity, iHeggie fe Rasiol (119961 ) 
provided an analytic estimate of this perturbation for a 
circular binary system in parabolic orbit around a third 
body. For the case of M. M b , the perturbation reads 
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with i the inclination of the CO binary relative to the or¬ 
bit around the hole, and <p depending on the initial phase 
of the binary and the longitude of the ascending node. 
Inserting the definition of tidal radius, Ea. Il2l becomes 


The following parameters in the three-body system are 
kept fixed: SMBH mass (M. = 10 6 _Mq), CO binary 
masses (mi = m 2 = m = 10 M 0 ), initial binary ec¬ 
centricity (eo = 0 ), and initial binary semi-major axis 
(ao = 10 Rq = 0.047 AU). With these parameters, the 
CO binary has a period 
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where M 6 = M./ 10 6 M 0 and r g = GM./c 2 is the gravi¬ 
tational radius. 

The parameters we vary are the CO binary orbital in¬ 
clination i, the initial longitude of the ascending node Ho, 
the binary initial phase do, and the pericentric distance r p 
via the penetration factor (3. Thus, our simulations of CO 
binary encounters with a SMBH span a four-dimensional 
parameter space {/3, Ho, do}- We sample this parameter 
space with random values taking uniform distributions in 
/3 _1 £ [0.35,5], cost £ [—1,1], and H 0 £ [0, 27t] . For the 
phase d 0 , we take 200 evenly spaced values between 0 
and 27 t. The distribution in /3 _1 implies that the values 
of the pericentric distance are uniformly distributed in 
rt £ [0.865, 8.65] AU. Furthermore, since the duration of 
the encounter is of order 
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the range of parameters we consider for /3 _1 imply that 
T p £ [0.1, 3] Pbp. Moreover, in the nomenclature of Heg- 
gie, our binary system is “soft” since 

\M b V. p 2 > -E bfi (17) 


which can be seen from 


Se = 6 y/n 2 1,/4 /3 3 ^ 4 exp 
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(13) 


This expression will be compared against our simulation 
results in section 16.11 Since Se vanishes as inclination 
t —X 7 r and /3 —> 0, we expect to see stronger agreement 
between simulation results and Eq. [13] for the low incli¬ 
nation and large pericenter encounters. 


3. BINARY ENCOUNTER SETUP 

The CO binary is injected in a parabolic orbit around 
the SMBH. The integration runs until either the binary 
is tidally disrupted by the SMBH and an amount of time 
equal to the initial time passes, or the center of mass 
of the binary has reached a true anomaly of O = — ©o, 
where @0 < 0 is the initial true anomaly of the orbit of 
the binary’s center of mass about the SMBH. 
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Also, since r t ~ 173r g , even for the deepest penetration 
encounters with /3 = 2, the pericentric distance will be 
r p = 87 r g . Therefore, it is perfectly safe to ignore gen¬ 
eral relativistic effects in all our encounters. Similarly, 
given that at pericentric passage 


w 2 GM. 

c 2 c 2 r p 
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post-Newtonian corrections to the orbital motion of the 
binary are at the level of a percent and thus will be ig¬ 
nored. However, in analyzing the the aftermath of an 
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Fig. 1. - K DE probability distribution function / as computed 
from Eg. l2ll for each of the {/3, l, #o} parameters. The pdf / 
represents the probability that a binary with a particular end-state 
began with a particular value of a certain parameter. The lines- 
class correspondence is: solid denotes class DB, dot-dash class long 
BEMRI, dash class short BEMRI, and fine dash class SU. 

encounter, we will take into account the Peters lifetime 
r gw of the binary, if survived, or if the EMRI if the bi¬ 
nary is disrupted. For the CO binary systems that we 
consider, the Peters lifetime du e to the emiss ion of GWs 
in the absence of the SMBH is ([Peters!1 1964) 


Tg W ,o = 0.95 x 10 8 yr 


a 0 


M b 

10 R@ ) 20 Mq 


-3 


( 20 ) 


The present study consists of N c « 13 million individ¬ 
ual simulations, resulting in: 2.1 million simulations of 
DB type encounters, 1.7 million yielding long BEMRI, 2 
million producing short BEMRIs and 7.1 million of the 
SU type. We integrate the Newtonian equations of mo¬ 
tion for the three-body system using Burlish-Stoer ex¬ 
trapolation with a leap-fro g integrator as described in 
iMikkola fe Tanikawal (|1999f ). We increase the accuracy 
of the i ntegration by the use of the CHAIN concept ala 
IMikkola fe Aarsethl (I1993T ). The coordinate system used 
for the numerical integration are barycentric coordinates 
with the three body center of mass at the origin. The 
orbital angular momentum of the CO binary-SMBH sys¬ 
tem is aligned with the z-axis. That is, the plane of the 
CO binary center of mass orbiting the SMBH is the xy- 
plane. We inject the CO binary in the first quadrant 
(x,y > 0 ) at a distance 200 r p and orient the orbit of 
the CO binary about SMBH such that pericentric dis¬ 
tance r p occurs along the ar-axis, specifically at y = 0 
and x < 0. Conservation of energy and angular momen¬ 
tum is checked at every time step, and simulations are 
halted and rejected if either quantity deviates from the 
initial value by one part in 10 . 


4. PARAMETER PROBABILITY DISTRIBUTIONS AND 
TIDAL RADIUS 

As already mentioned, we have classified the encoun¬ 
ters into four types: DBs, long and short BEMRIs and 
SUs. Next, for each of these outcome types, we estimate 
the probability distribution as a function of the varying 


parameters {/?, i, D 0 , $o}- For example, f(p~ 1 \DB) is the 
probability distribution function for the parameter p~ x 
for disrupted binaries, i.e. the probability that a dis¬ 
rupted binary had a particular value of /3 —1 (solid curve 
in upper left panel of Figure [1]). We use for this purpose 
the technique of kernel density estimation (KDE). KDE 
is a non-parametric method for estimating probability 
densities in which a kernel function K is convolved with a 
collection of Dirac delta functions. KDE asymptotically 
converges to the true distribution faster than histograms 
dScottl I1979h . The KDE probability distribution / of a 
parameter x is computed from 

1 N , N 

f 0) = ^2K(x)*5(x-x n ) = — K{x-x n ) (21) 

n —0 n —0 

where N is the number of data points in the sample. We 
use a Gaussian kernel K{x) = (27rh 2 ) _1//2 exp [a; 2 /(2/i 2 )] 
with a variance of h 2 = [parameter range]/100, chosen 
to produce distributions that retain structure while not 
being over-smoothed. Figure [T] shows the resulting nor¬ 
malized KDE probability distributions / for each of the 
{/?, t, D 0 , $o} parameters obtained from Eq. [21]. 

We can draw from these probability distributions sev¬ 
eral conclusions about the nature of the end-state of the 
binary relative to the input parameters, as well as the 
predictive power of the individual parameters. It is clear 
from the panel of the p^ 1 probability distribution (top 
left panel in Figured]) that the encounters yielding binary 
disruption (i.e. DB-type, black, dot-dash line) occur pri¬ 
marily for small values of /3 _1 or large penetration factors 
/ 3 . This should be obvious: closer passes translates into 
stronger tidal forces and higher likelihood of disruption. 
Additionally, regarding the parameter l, prograde bina¬ 
ries (defined here as those with inclination t < n/2) are 
more probable to disrupt than retrograde binaries (those 
with i > 7 t/ 2 ). This is so because, in general, retrograde 
binaries are more likely to survive and become bound to 
the SMBH, implying that while the BH orbit loses en¬ 
ergy, the surviving retrograde binaries must gain an equal 
amount energy, though not generally enough to disrupt. 

Notice that the probability distribution for the initial 
orbital phase 9q is flat. That is, the value of 9q does 
not play a role on the outcome of the encounter. The 
distributions in Figure [21 are over all the simulations in 
the parameter space. The flat distribution in 9q illus¬ 
trates that there is no strong correlation between 9q and 
disruption across the entire parameter space of interest; 
the combination of other parameters has far greater in¬ 
fluence. Clearly if all parameters except 9q were fixed, 
some binaries would be disrupted, but there is no con¬ 
sistent trend. Also interesting is that the probability 
of longitude of the ascending node Do is relatively flat. 
Therefore, we will focus our attention on the parameters 
P and l since they have the largest effect on the binary 
end state. 

The fact that p and l are the most relevant parameters 
can also be seen in Figured] where we show the fraction 
of surviving binaries for each parameter. Consistent with 
FigureQ] the top panels in Figure [2] show that the sur¬ 
vival probability depends more strongly on /3 _1 and i. 
The top right panel shows again that prograde (i < tt/2) 
binaries in general more likely to be disrupted than ret- 
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Fig. 2.— Fraction of surviving binaries as a function of each in 
parameter for each parameter {/}, Qq , Oq }. 


rograde (t > 7 t/2 ) binaries. Similarly, it is clear from the 
top left panel that the probability of disruption goes to 
zero for /3 _1 > 2.1. Also consistent with Figure [U the 
bottom panels in Figure [5] show that fl 0 and 9 0 have a 
very weak influence on the survival ratio. Therefore, we 
can assume with confidence that for practical purposes 
the parameter space is two-dimensional, namely 
The dependence of the binary survival on i implies that 
the tidal radius does not only depend on f3 but also on the 
inclination of the binary. This is not surprising because 
the forces from the SMBH responsible for disrupting the 
binary are those projected along the semi-major axis of 
the binary. Figure [3] depicts a two-dimensional histogram 
of tidal disruptions as a function of cost and /3 _1 . The 
solid line is the condition obtained by setting in Eq. [12] 
the eccentricity perturbation to Se = 1 with <j> = 0 to 
maximize the effect. Namely, 


.F - 1 (<,, 0 ) = 6 -v/ 7 r 2 1 / 4 / 9 _3/4 exp 


fL^o-3/2 


( 22 ) 


Notice that this analytic approximation bounds the re¬ 
gion containing disruptions well for /3 _1 > 1 , but does 
not do so well for < 1 , where the deepest penetration 
encounters are located. 


5. DISRUPTED BINARIES 

As stated in the introduction, CO binaries disrupted 
by the SMBH are of interest as both sources of EMRIs 
and HVSs. In this section, we investigate the channels 
for EMRI formation and use previous HVS results as a 
check for our simulation accuracy. 

We recall that after a binary is disrupted, the total 
energy of the system is approximately given by E = E\ + 
E% = Ebfi, where E ^0 is the initial binding energy of the 
CO binary and 


1 2 GrriiM, 

Ei = 2 miVi --- 


(23) 


with i = 1 , 2 . 

Figure [4] shows the combined histograms of the normal¬ 
ized energies Ei/\Eb^\ and -E^/I-E^ol (top panel) and the 
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Fig. 3.— 2D histogram of disruptions vs. /3 1 and cos l, sho wn 
as fraction of total disruptions. The white line denotes Ea. 1221 
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F ig. 4.— Top panel: Histogram of the binding energy Ei from 
Eq. 1251 normalized to the initial binding energy of the binary E 5 q- 
Middle panel: Histogram of the corresponding eccentricity. Bottom 
Panel: Histogram of resulting EMRI semi-major axis. Histograms 
shown as fraction of total disruptions. 


corresponding eccentricities (bottom panel). The sym¬ 
metry of the histograms gives the impression when a bi¬ 
nary is disrupted the outcome is invariable one compo¬ 
nent bound to the SMBH, i.e. an EMRI, and the other 
component unbound to the hole, namely a HVS. This is 
indeed the case for most of the disruptions (~ 99.9967%). 
However, since E\ + E 2 = E^, 0 < 0, there are cases, al¬ 
though rare (~ 0.0033%), for which both E\ < 0 and 
E 2 < 0, and the outcome is two EMRIs. Moreover, no¬ 
tice from the bottom panel in Figure [3] that the average 
eccentricity for the bound (unbound) component of the 
disrupted binary is e_ ^ 0.97 (e + ^ 1.3). This can be 
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understood from the definition of eccentricity 

2 E±L 2 ± 


6± 1 ' G 2 mi M 2 ’ 


where we approximate 


E± ft ±G M. to± a 0 r 


and 

Thus, 
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Furthermore, from Eg. E5l and E b}0 = |GM 2 /(4a 0 )|, we 
have that 


E± 


I E, 


6,01 


±74/3 2 M 6 1/3 


m 6 

20 M 0 


-1/3 


(28) 


consistent with the histograms in the top panel of Fig¬ 
ure m 


5.1. Ejected Hyper-velocity Stars 

Hypervelocity stars are stars with velocities of the or¬ 
der of hundreds or thousands of km s _1 , which may 
exceed the escape veloc i ty of our galaxy. These stars 
were predicted by iHilld (l'l988t l as the result of binary 
disruption in the galactic center, an d discovered obser- 
vationally nearly two decades later dBrown et al.l [20051 : 
lEdelmann et al.ll2005f) . 

The production of HVS from the tidal disruption of 
binary systems by a SMBH ha s been discussed b vlHilkJ 
(11988h: I Y u fc Tr emainel (|2003H : lAntonini et al.l (|2010H : 
IBromlev et al.1 (I2006T) . Our study is consistent with those 
results. From Eq. l25land Ei = rm 4/2, one has that 


4 ~ 2 GM.a 0 r p 2 
ft 2 G /3 2 M. a 0 r t “ 2 

ft 2 G /3 2 Ml /3 a/" 1 M 2/3 , (29) 


which yields the following asymptotic velocity approxi¬ 
mation of a HVS 


i 5, 312 km/s -/3 -M 6 


1/6 


(a 


a 0 


047 AU 7 


n-i/2 / M b \ 1/3 


20 M 0 


(30) 


Resc al ed to our case, the numerical study by 
IBromlev et al.i (|200fil ) found that 


Voo ~ 4,468 km/s • g(D) ■ M^ 6 

x ( a ° Y 1/2 ( Mb V /3 

X V0.047AU7 \20M q ) 

where 


g(D) = 0.774 ± 0.0204 D - 6.23 x 10~ 4 D 2 
± 7.62 x 10~ 6 D 3 - 4.24 x 10“ 8 D 4 
±8.62 x 10 ~ n D 5 . 


(31) 


( 32 ) 




(1000 km/s) 


Fig. 5. — Top panel: Grey dots denoted the ejection velocity 
Voo as a function of f3 for all the unbound stars from disrupted 
binaries. Stars show the averages over 100 evenly spa ced intervals 
in (3 . Solid line shows the model prediction from Eq. l3ll Bottom 
panel: distribution of i?oo for all ejected components. 

with D = 79.37 /3 _1 . 

In the top panel in Figurc[5]we show with grey dots the 
ejection velocity Voo as a function of /3~ 1 for all the un¬ 
bound stars from disrupted binaries. In the same panel, 
stars show the averages over 100/3 bins. The solid line 
shows the model prediction from Eq. 1311 T he agreement 
of our average velocities with those from IBromlev et al.1 
(f200l is evident. However, it should be noted that the 
range of possible velocities can vary significantly from 
this average, with the largest values roughly double the 
analytic prediction and nearly 59% of ejections in our 
simulations exceeding the predicted value. 


5.2. Extreme-Mass-Ratio Inspirals 

The traditional formation channel of an EMRI is when 
two-body relaxation not only bounds a star to a SMBH 
but also brin gs i t into the G W-emission inspiriting regime 
(lAmaro-Scoane et al.l 12007 11. The pericentric distance to 
capture a star bound to a SMBH with semi-major axis 
a c and form an EMRI is 


3 r g M 6 


-4/7 


TO 


10 Mr. 


2/7 


0.05 pc 


2/7 


(33) 


Since the star needs to penetrate so deeply into the the 
neighborhood of the SMBH, despite circularization from 
GW emission, the EMRI arrives to merger with signifi¬ 
cant eccentricity, e ~ 0.5 — 0.9. 

The disruption of a binary by a SMBH provides an al¬ 
ternative channel since the disruption will always leave 
at least one component bound to the SMBH. The cap¬ 
ture distance in this situation is much larger. It is not 
r c but the tidal radius r t ~ 173 r g (see Eq. [15]). Being 
bound to the SMBH does not necessarily translates into 
becoming an EMRI. The condition for a star to qual¬ 
ify as an EMRI is that the timescale r gw for its orbital 
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Fig. 6 .— 2D histogram of potential EMRI candidates as a function of eccentricity and semi-major axis a^, shown as fraction of total 
disruptions. White lines are lines of constant Peters lifetime with values r gw = 10 a years, where from left to right a = 5.0, 5.5, 6.0, 6.5, 
7.0, 7.5, 8.0 and 8.5, respectively. Red lines show lines of constant /3 ~ 1 with /3 ~ 1 = 0.35, 0.5, 0.75, 1.0, 1.25 1.5 1.75, and 2.0 from left to 
right. 


decay by GW emission is sufficiently shorter than the 
relaxation t i me £ r ix of the stars in the galactic center 
(lAmaro-Seoane et al.ll2007h : 


1-8 x 10 s yr 


10 Mq 


( _z_V 

VlOOkms 1 / 
lO 6 M 0 pc -3 


(34) 


will have a typical eccentricities of eo ~ 0.97 and semi¬ 
major axis cio ~ 45 AU; thus, cq ~ 2.3 AU in Eq.[55] For 
a space-based interferometer like LISA, the low end of 
the sensitivity window is / gw ~ 10~ 4 Hz. If we approx¬ 
imate / gw « 2/ with / _1 = P = 2 7r a?/ 2 /y/G M . the 
Keplerian orbital period, we have that when the EMRI 
enters the sensitivity band of a LISA-like detector it has 
a semi-major axis of 


where cr is the local velocity dispersion, n is the local 
number density of stars, and m is the average stellar 
mass. 

Figure [6] displays a 2D histogram of (a, e) for EMRI 
candidates from the disruption of binaries. White lines 
are lines of constant Peters lifetime with values r gw = 10“ 
years, where from left to right a = 5.0, 5.5, 6.0, 6.5, 7.0, 
7.5, 8.0 and 8.5, respectively. Red lines show lines of 
constant /3 _1 with /I -1 = 0.35, 0.5, 0.75, 1.0, 1.25 1.5 
1.75, and 2.0 from left to right. It is clear from this 
histogram that most of the candidates qualify as EMRIs 
since r gw -C t r i x ps 10 8 yr. Also interesting is that the 
encounters with deepest penetration factors, i.e. small 
/3~ 4 , yield EMRIs with larger eccentricity. 

Our results are consistent with those from lMillcr et al.l 
(l2005f) in the following respect. Because their “capture” 
radius is much larger than in the traditional scenario, 
i.e. rt ~ 170 r g instead of r c ~ fewr ff , the EMRI from 
binary disruptions will have circularized dramatically by 
the time they merger or they enter the sensitiv e band 
of a L ISA-like detector. One can see this using I Pc ter 4 
<1961 


a(e) 


co e 12 / 19 / 121 2 \ 870/2299 

(1 -e 2 ) V + 304 6 ) 


(35) 


with Co determined by the initial condition a = ao when 
e = eo- From Figure^] we have that at “birth” the EMRI 


a < 74 r 0 




gw 


10- 4 Hz 


-2/3 


M, 


-2/3 


(36) 


From Eq. 1351 one has that the corresponding eccentricity 
is e < 0.15. By the time the EMRI merges, i.e. a ~ r g , 
it will have eccentricity e ~ 10 -4 and emit GWs at dom¬ 
inant frequencies of / gw ~ 60 mHz. This is clear from 
Figure [T] where we plot the eccentricity of the EMRIs 
as a function of the GW frequency / gw . We approxi¬ 
mate /gw « 2/ with / _1 = P = 2 7r a 3 / 2 A /G AI, the 
Keplerian orbital period and a from (lPeterslll964f ). Fig¬ 
ure [7] shows that in the sensitive LISA band (0.1 mHz < 
fGW < 100 mHz), EMRIs can have average eccentricity 
e < 0.1. At the most sensitive frequencies, however, at 
few 10 mHz, these EMRIs will have typica l eccen tric- 
ities e < 0.01, in agreement with lMiller et al.l (|2005f l. 


6 . SURVIVING BINARIES 

We now turn the attention to the population of bina¬ 
ries that survive the encounter with the SMBH. These 
binaries make up the majority of the encounters we sim¬ 
ulated given the range of values of j3 we considered. In 
these scenarios, the energy exchange is simpler: energy 
is either drained from the orbit of the binary around the 
BH and donated to the CO binary orbit, or vice versa 
(see Eq. [6]). In general terms, we expect the binary to 
soften and become bound to the SMBH or to tighten and 
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^ gw 


Fig. 7. — EMRI eccentricity at fundamental frequency / gw . Ec¬ 
centricity for individual EMRIs are shown in gray, average shown 
in black. 

remain unbound. Specifically, the surviving binaries will 
belong to one of the following types: binaries bound to 
the SMBH with a Peters lifetime r gw longer that their 
orbital period P, about the hole (i.e. long BEMRIs); 
binaries also bound to the SMBH but that merge be¬ 
fore completing an orbit around the hole, r gw < P . (i.e. 
short BEMRIs); and binaries that are unbound from the 
SMBH. We will analyse the distributions of orbital pa¬ 
rameters (eccentricity and semi-major axis) of the surviv¬ 
ing binaries, with the goal of determining the net change 
to merger lifetime and whether this leads to detectable 
changes in the CBC rate observed by detectors such as 
LIGO. But before doing so, we will briefly revisit the 
effect of and Ho- 

In Section U we estimated probability distributions of 
the parameters {/3, l,9q,^o} for each of the encounter 
types. We observed that the distributions for the phase 
of the binary do and the longitude of the ascending node 
Ho are to a good approximation flat. Thus, the param¬ 
eters (3 and l have a more dominant role on the out¬ 
come. This does not mean that choices of do and Ho 
do not influence the end state. Figure [8] displays the 
eccentricity of the CO binary after the encounter with 
the SMBH as a function of do- In the upper panel, 
we show with dots cases with parameters (/3 —1 , t. Ho) = 
(0.77,123°, 46°), squares with (1.9,101°, 312°), and di¬ 
amonds with (4.6,49°, 152°). Points where the curve 
jumps above 1 indicate disrupted binaries. Notice as 
expected that as the penetration factor increase, so the 
eccentricity gained by the binary as well as its depen¬ 
dence of the initial binary phase. In the lower panel in 
Figure [8] we show with circles parameters (/3 —1 , 6, Hq) = 
(0.57,177°, 255°) and with squares have (1.7, 79°, 308°). 
Notice that binaries making a closer pass (circles) gener¬ 
ally suffers a reduced perturbation to eccentricity due to 
the nearly retrograde orbit. 

6 .1. Post-Encounter Binary Eccentricity 



Fig. 8. — Resulting eccentricity of surviving binaries as a func¬ 
tion of Oq. Upper panel: the dots show cases with parameters 
(/3~ = (0.77,123°, 46°), squares with (1.9,101°, 312°), and 

diamonds with (4.6,49°, 152°). Points where the curve jumps 
above 1 indicate disrupted binaries. Lower panel: dots with 
parameters (/3 — 1 ,l, f2o) — (0.57,177°, 255°), and squares with 
(1.7, 79°, 308°). 
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Fig. 9. — Fraction of surviving binaries as a function of post¬ 
encounter eccentricity as well as with a solid line the cumulative 
distribution function (CDF) of the eccentricity values. 


In Figure [5] we show the fraction of surviving bina¬ 
ries as a function of post-encounter eccentricity as well 
as the cumulative distribution function (CDF) of the ec¬ 
centricity values. Note that the majority of the surviving 
binaries are relatively unperturbed in eccentricity, with 
the 65% quantile lying at approximately e ~ 0.109 and 
the 90% quantile at e & 0.553. 

Figure [TO] show how the post-encounter eccentricity 
depends on /3 _1 and l. The surviving binary acquires 
eccentricity for larger penetration factors (smaller /3 _1 ) 
and prograde orbits (cos^ > 0). Contours of constant 
loge = —5.5, —5.0, —4.5,..., 0 from right to left are shown 
in white. Black lines are the corresponding estimates ob¬ 
tained from Eq. [13] averaged over (j). Notice that the 
contour for loge = 0 is not present for the simulation 
data. All simulation values have been averaged over 9q. 
As expected, we see better agreement between the pre¬ 
diction from Eq. Il3l and simulation results for large /3 _1 
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Fig. 10. — Resulting eccentricity loge for surviving binaries plot¬ 
ted as an intensity map vs. (3~ 1 and t. Contours of constant 
loge = —5.5, —5.0, —4.5, ...,0 are show n fro m right to left in white 
and for the analytic estimate from Eq. ECU in black. 

and small l since that is the regime where Eq. [13] was 
derived. The small overestimation of the prediction at 
small inclination is co nsistent with what was found in 
iHeggie fe Rasiol (Il996f) for an equal mass binary. 

Finally, Figure fill shows gray scale maps of the post¬ 
encounter loge against pairs of parameters /3 _ 1 ,t, and 
n 0 , averaged over 6q. Noticeable are the different effects 
of the input parameters. As /3 _1 increases, the final ec¬ 
centricity decreases toward zero, however even binaries 
with /3 _1 = 3 can receive a non-negligible amount of ec¬ 
centricity, e ~ 0.2, after the encounter. Inclination plays 
a smaller but important role in determining the allowable 
range of final eccentricities, with eccentricity decreasing 
as the binary orbit approaches retrograde with respect 
to the SMBH orbit. As stated before, the longitude of 
the ascending node flo plays no predictable role in the 
perturbed eccentricity. 

6 .2. Post-Encounter Binary Semi-major Axis 

Semi-major axis of the surviving binaries is shown as 
a histogram in Figure [TUI The semi-major axis has been 
normalized to its initial value. It is clear that the semi¬ 
major axis is generally not perturbed to a strong degree, 
with the overall mean being (a/ao) = 1.03. If the sample 
size is restricted to binaries with final semi-major axes 
in the range 0.5 < a/a o < 1.5, which accounts for > 97% 
of the data and removes the outliers, then the mean is 
(a/ao) = 0.988 and standard deviation of 0.084. 

Figurefl3lis similar to Figure fill but for a/ao. It can be 
seen here that for /3 _1 < 2 the binary can suffer a large 
change to the semi-major axis, while above this value 
the change is small. Inclination has the general effect of 
tightening the binary for i < 7r/2, and loosening it for 
l > 7t/2. There is some variation in the extent of semi¬ 
major axis perturbation as a function of flo, particularly 
when the inclination is around l ~ 7r/2 or at low values 
of /3 , with minimum change at flo « 0 or 7r. 

Since the binding energy of the binary is related to the 
semi-major as Eb oc — a -1 , one can from Figure [13] iden¬ 
tify regions of parameter space where the binding energy 
of the binary increases or decreases. Furthermore, given 


that E = E cm + Eb = Eb t o < 0, an increase in binary 
energy Eb must be compensated by a decrease in the 
BH orbit energy E cmj resulting in capture of the binary. 
Therefore the points in Figure [13] where a/ao > 1 corre¬ 
spond to binaries which become bound, while a/ao < 1 
indicates binaries that remain unbound. 


7. PETERS LIFETIME AND CBC RATES 


As mentioned before, the Peters lifetime is an estimate 
of the time that a binary system with eccentricity e and 
semi-major axis a t akes to merge as it looses energy due 
to emission of GWs (|Peterslll964f) . As presented in Equa¬ 
tion [20] the CO binary system in our study, with initial 
semi-major-axis ao = 10 Rq and vanishing eccentricity 
(i.e. circular orbit), has a Peters lifetime of r gWi0 « 10 8 
years. Figure [TJ] shows the Peters lifetime r gw normal¬ 
ized to T gWi o for all of the surviving binaries. Notice the 
similarity of the maps with those for the semi-major axis 
in Figure [12] 

Unbound binaries that survive the encounter with the 
SMBH, what we call SU binaries, could in principle have 
a shorter ( T gw < T gWj o) or longer (r gw > r gWi o) Peters life¬ 
time. The last situation, however, is not possible given 
our initial setup. An unbound binary requires increas¬ 
ing the energy of its center of mass E cm . This will come 
at the expense of decreasing its binding energy Eb since 
E = E cm + Eb = Ebfi < 0, which in turn requires de¬ 
creasing its semi-major axis. A decrease in its semi-major 
axis translates into a shorter Peters lifetime. Thus, all of 
the unbound binaries in our study will have accelerated 
merger times due to the SMBH 

We recall that a short EMRI is one for which r gw < P.. 
That is, the binary is expected to merge before returning 
to the hole. On the other hand, a long BEMRI is one in 
which r gw > P„. These BEMRIs have typically highly 
elliptical orbits with average eccentricity e cm ~ 1 — 10~ 5 
and long periods. Therefore, 


o 3/2 

2 7 T a C m 


2 7r r 


3/2 


P- ‘ , ',“C°1 - ~'"P 3/2 

Vgm;- VGW. ’ 


=r 3 / 2 (i-e cm )- 3 / 2 2 7rr ‘ 


3/2 


VGW. 


= r V \ l-e c 


o 3/2 
\— 3/2 2 7 T a 0 
v 


= 7.2 x 10 4 yr/3 " 3/2 

-3/2 

/ J- C-rm \ 

X 


10" 5 


VGW b 


Pb, 0 

7.16 x 10 4 s 


(37) 


The CO binary in a long BEMRI could potentially sur¬ 
vive several orbits around the SMBH, with its orbit being 
perturbed by each passage. We will investigate the fate 
of BEMRI systems will be the focus of a future study. 

Independently of being in an BEMRI or unbound, CO 
binaries that after an encounter with a SMBH are af¬ 
fected in such a way that r gw < r gW! o have particular in¬ 
terest. They could potentially increase the rates of CBC 
to be detected by LIGO. We construct a simple formula 
to obtain a first estimate of this effect as 


Scbc = [r * fb* Nc(D h )] * ( St * / l ), (38) 

where Scbc is the enhancement to the current predicted 
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Fig. 11.— Gray scale maps of post-encounter loge as a function of combinations of input parameters. All values have been averaged over 
00 - 



Fig. 12. — Histogram of the binary semi-major axis distribution 
after SMBH encounter. The horizontal axis has been stretched 
so that the few outliers did not dominate the plot scale, and the 
vertical axis has been truncated to show detail with the central 
bin having a true value of ~ 0.69. It is clear that the majority of 
encounter result in very small change in semi-major axis. 

CBC rate Rcbc , be., Rcbc = Rcbc + £cbc • The 
first three terms in Eq. [38] are observational quantities 
taken from the literature. T is the estimated encounter 
rate of stellar mass objects / per SMBH per year, /& is 
the binary fraction, and Nc{Dh) is the number of Milky 
Way equivalent galaxies (MWEGs) observable by a GW 
detector with horizon distance, Dh. The product of these 
factors is the rate of binary encounters with a SMBH in 
the observable volume of a GW detector. 

The next two factors come from our simulation 
results. The single binary enhancement £ T = (1 — t/to) 
is the percent difference between the old merger 
lifetime, To, and the mean new lifetime, f, and 
fz is the fraction of binaries from our simulations 
that result in guaranteed LIGO sources after the 
SMBH encounter. Encounter rates in galaxies with a 
1O 6 M 0 SMBH, estimates for the single CO capture 
rate range from T ~ 5 x 10 -9 yr -1 MWEG -1 to as 
high as T ~ 10~ 6 yr" 1 MWEG -1 dHils fc Benderl 

(119951): iSieurdsson fc Reei (1199711: llvanovl (120021) : 
iHonman fc Alexander! (120051) : iMerritt et al.1 (120111) ) ; 
we use the higher side of these estimates, 
T = 10 -7 yr -1 MWEG -1 , as the pericenter dis¬ 

tances in our simulations reach much larger values than 
the single star capture radius. The binary fraction, fb, 
near the galactic center is not a well determined value, 
however, in the absence of better knowledge, we take it 
to be roughly the same as the binary fraction of field 


stars, fb = 0.5. The number of MWEGs observable 
by aLIGO with a BH-BH mer ger horizon distanc e 
of Dh = 2187 Mpc is given as (|Aba,di e et al.1 (|2010f) : 
iKaloeera et al.l (120011) ) 


NG = t v 

G 3 V M PC/ (2.26) 3: 


(39) 


which gives a value of Ng ~ 4.4 x 10 7 MWEGs. 

From our simulation results, we find that binaries from 
the fz, categories have a mean Peters lifetime of f « 
0.84r o giving £t ~ 0.16 and LIGO fraction fz « 0.68. 
Combining these factors, we compute an estimated en¬ 
hancement to the CBC rate of £cbc ~ 0-25 yr -1 . The 
predicted rate of expected BH-BH mergers has been es¬ 
timated to lie between 0.4 MWEG -1 Myr -1 < T BH < 
30 MWEG -1 Mvr -1 f or realistic to o ptimi stic scenarios 
iKalogera et al.l ( 20071) : lAbadie et all (|2010l ). This corre¬ 
sponds to an estimated merger rate within the aLIGO 
volume of ~ 20 yr -1 < N B h < 1300 yr -1 , which our es¬ 
timated rate enhancement £cbc may change by as much 
as ss 1%. This enhancement will likely be difficult to de¬ 
tect with small observation catalogs and given the uncer¬ 
tainty in the estimated merger rates, but could become 
noticeable over long observation times. 

8 . SUMMARY AND CONCLUSIONS 


In this paper, we have presented the results of ~ 13 
million individual simulations of parabolic encounters be¬ 
tween a CO binary and a galactic center SMBH while 
varying the orientation of the binary and its distance of 
closest approach to the SMBH. 

Tidal disruption of the binary occurs with about 16% 
probability given the full range of parameters covered in 
this study. Consistent with previous work in this area, 
this sort of disruption can create HVS which can escape 
from the SMBH with high speed. We also explored dis¬ 
ruption as a formation mechanism for EMRIs, which are 
of interest to space-based GW detection missions, and 
found that the EMRIs formed in this way will gener¬ 
ally have very low eccentricity when they enter the LISA 
band. This work shows that considering the full range of 
possible orientations gives a broader range of formation 
eccentricities than previous estimates have predicted. 

Surviving binaries can either become bound to the 
SMBH after the encounter and become a BEMRI or re¬ 
main unbound from the hole. Among them, there is a 
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Fig. 13.— Resulting semi-major axis (a/ao) for surviving binaries plotted as intensity maps vs. pairs of input parameters. All values 
have been averaged over Oq. 
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Fig. 14.— Gray scale maps of the Peters lifetime log(rg W /Tg W ,o) for surviving binaries vs. pairs of input parameters. All values have 
been averaged over 6q. 


subclass of surviving binaries for which r gw < r gWj o those 
with a small enough BH orbital period such that they will 
not merge before completing one orbit are the “short pe¬ 
riod” BEMRIs, and understanding their full evolution 
requires a more careful (i.e. Post-Newtonian) approach 
to the integration in order to account for eccentricity and 
semi-major axis change due to GW loss during the long 
orbit. Both the unbound binaries and the long period 
BEMRIs which will merge before one SMBH orbit are 
factored into the calculation for the CBC rate enhance¬ 
ment, which is potentially important to ground based 


GW detectors such as LIGO. We find that for the aLIGO 
volume, the enhancement factor for BH-BH mergers is 
£cbc ~ 0.25 yr -1 or as much as 1% of the predicted 
rates, though this enhancement may be difficult to sta¬ 
tistically detect in accumulated event catalogs. 
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